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Rocket Based Combined Cycle (RBCC) engines operate within a wide range of Mach numbers and altitudes. 
Fundamental fluid dynamic mechanisms involve complex choking, mass entrainment, stream mixing and wall 
interactions. The Propulsion Research Center at the University of Alabama in Huntsville is involved in an on- 
going experimental and numerical modeling study of non-axisymmetric ejector-based combined cycle 
propulsion systems. This paper attempts to address the modeling issues related to mixing, shear layer/wall 
interaction in a supersonic Strutjet/ejector flow field. Reynolds Averaged Navier-Stokes (RANS) solutions 
incorporating turbulence models are sought and compared to experimental measurements to characterize 
detailed flow dynamics. The effect of compressibility on fluids mixing and wall interactions were investigated 
using an existing CFD methodology. The compressibility correction to conventional incompressible two- 
equation models is found to be necessary for the supersonic mixing aspect of the ejector flows based on 2-D 
simulation results. 3-D strut-base flows involving flow separations were also investigated. 


I. INTRODUCTION 


^^fundamental understanding of ejector flow physics is an enabling technology to realize an operational Rocket 
^Based Combined Cycle (RBCC) propulsion system. Past theoretical and experimental ejector studies have 
considered one-dimensional, axisymmetric, or at best two-dimensional geometries. Concepts such as the Strutjet 
use a complex asymmetric, three-dimensional flow path. The Propulsion Research Center (PRC) of University of 
Alabama in Huntsville (UAH) has an ongoing research program to characterize asymmetric ejector performance in 
terms of mass flow entrainment, stream mixing, and choking mechanisms using both an experimental approach and 
computational fluid dynamics (CFD) modeling approach 1 

Validation of computational fluid dynamic (CFD) models appropriate for supersonic through hypersonic flow 
applications requires careful consideration of the physical processes encountered at these flight regimes and detailed 
comparisons of the calculated results with experimental data sets that include these processes. Supersonic ejectors 
perform the mixing and recompression of the primary flow (ejector) and the secondary (induced) flow 2 . The primary 
flow gets accelerated by a convergent divergent nozzle to reach supersonic velocity. By an entrainment-induced 
effect, the secondary stream is drawn into the flow and accelerated. Mixing and recompression of the resulting 
stream then occurs in a mixing chamber, where complex interactions take place between the mixing layer and 
shocks. In Reference 1, several turbulence models ranging from the algebraic to two-equation models were 
systematically evaluated for 2-D single-nozzle ejector flow field. Depending on the operating conditions of the 
ejector chamber pressure, the mixing patterns and the shear layer/wall interactions exhibit drastic changes. In 
general, it was concluded that the SST two-equation model of Menter gave the most satisfactory results 1 . The SST 
model is a combination of two popular two-equation models: the k - co and the k-z models. A blending function was 
devised such that k - co model was used in the wall region and a switch to the k-z model occurred in the outer 
wake/ffee stream region. In supersonic mixing regime, the k-z model was found to over-predict the mixing of 
compressible shear layers 3 . Various compressibility corrections to the k-z model was proposed [see Chapter 5 of 
Reference 4], and the performance of these models were inconclusive. 
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For the flow field encountered in the present study, it was not clear what role compressibility corrections would 
play. The purpose of this paper is to incorporate an existing compressibility correction into the k-z model, for 
calculating mixing patterns of ejector flows under various operating conditions. In addition, the 3-D features of the 
investigated flow will be discussed and modeled. 

II. NUMERICAL MODELING 


Compressibility effects, which reduce mixing and turbulent intensity are present in supersonic free shear aero- 
propulsive flows and can be significant for large convective Mach numbers. A parameter that could be used to 
characterize the effect of compressibility on the turbulent fluctuations is the turbulent Mach number, which is 
defined as the ratio of a characteristic speed of turbulence to the wave speed. Thus M t = V2k/a, where ‘k’ is the 
turbulent kinetic energy and 4 a’ is the speed of sound. As pointed by various researchers 5 6 the standard k- 8 model 
fails to predict the observed decrease in spreading rate for the compressible free shear layers. Compressibility is 
modeled using dilatational and pressure fluctuation terms within the turbulence quantity conservation equations. 
Concentrating on the k- e model, the dissipation of the kinetic energy in compressible turbulence was argued to be 
augmented by the dilatational components 8d of the dissipation tensor in addition to the solenoidal components. 
Sarkar 7 et al., considered the evolution of the turbulence fluctuations on an acoustic time scale. They considered the 
effect of varying compressibility based on the turbulent Mach number M t , on the rate of dissipation of turbulent 
kinetic energy, e. Analysis of decaying compressibility on the solenoidal component of the dissipation rate, e s , was 
negligible in comparison with the effect on the dilatational component, 8d. They also observed that the ratio of 
dilatational dissipation rate to the solenoidal dissipation rate varied directly as die square of the turbulent Mach 
number M t . They proposed a model for the dilatational dissipation rate, which is given as 

Ed = (Xi £s M t 2 , where ai is an arbitrary constant of order 1. 


Adding the above term would account for the compressibility correction in the conventional two-equation model. 
The baseline two-equation model used in this study is the standard k-z model, written for mass weighted turbulence 
quantity for turbulence kinetic energy with the extra compressibility term being: 
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Sarkar et al. 7 and Zeman 5 have postulated that the transport equation for 8 is unaffected by compressibility. The 
compressibility corrections were incorporated into the density-based CFD-FASTRAN® code 8 , for the k-z turbulence 
model. 


ffl. RESULTS AND DISCUSSION 

Before modeling the supersonic ejector flows, the implemented compressibility model was first validated with 
the supersonic backward-facing benchmark case of Eklund et al 9 . It should be stated that grid-refinement studies 
were carried out 10 and the results presented here are grid- independent. 

Supersonic Backward Facing Step 

The Mach 2 flow over a backward-facing step has been investigated experimentally using a variety of 
techniques. The test section geometry and the reattachment length calculations are shown in Table 1. The flow field 
contains many complex flow features including boundary layer separation, reattachment, viscous shear layer and 
expansion fan/shock wave interaction. In the calculation, the expansion fan, centered on the step, turns the flow over 
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the step. Behind the step a recirculating region of low-momentum fluid that is characterized by low velocities and 
low pressures is formed. Immediately behind the step a shear layer forms between the accelerating flow turning over 
the step and the recirculating flow. Downstream of the step a reattachment shock turns the flow parallel to the wall. 

The net effect of the compressibility correction is to reduce the magnitude of turbulence kinetic energy, thus the 
effective eddy viscosity, and to delay the flow reattachment. Both the standard k-e and the compressibility corrected 
k-e model utilized the wall-function approach in this study. The SST model of Menter utilized £-co blending for near 
wall treatment and recovered the standard k-e model at the outer region. Thus the SST model is viewed as the wall- 
resolving k-e model without compressibility correction. Figures la and lb show the calculated velocity profiles 
(normalized with respect to mean inlet velocity) at the upstream and downstream locations of the reattachment. As 
can be seen from the figures, due to the reduction of shear layer mixing, the compressibility correction gave longer 


Table 1. Test Geometry & Reattachment Lengths 


Feature 

Dimension 


Reattachment Length 

Step Height 

H=0.003 18m 

Experimental 

3.9H 

Test Section Height 

6.71H 

Menter SST 

3.20 

Test Section Length 

10H 

Standard k-e 

3.65 



k-e /compress. 

3.81 


reattachment length and compared better with experimental data. The SST model, though it was originally based on 
k - co compressible flow formulation, did not have explicit compressibility correction at shear layer region and under- 
predicts the reattachment length. 
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Figure lb. Normalized velocity profiles at X=61.7H (downstream region) 
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UAH’s Strutjet Ejector 


As described in Lineberry et al., 11 the single-nozzle asymmetric ejector flow system involves complicated 
phenomena including shock-boundary layer interaction, mass choking, jet mixing and secondary flow separations. 
In Reference 1, the ejector mass choking phenomena created by a high-pressure jet within confined environment, as 
well as the mixing patterns were studied. The compressibility effect of interest here is applied to the mixing region 
(see Figure 2). Due to the large aspect ratio geometry of the single-strut configuration, it was beneficial to carry out 
detailed 2-D simulations (see also Reference 1) to first investigate the relative performance of the turbulence model 
with compressibility correction. The sidewall database was chosen for this purpose. The ejector geometry, grid 
systems, as well as the boundary condition specifications were described in Reference 1. 



Figures 3 to 5 show the comparison of calculated side-wall pressures with the experimental data for three 
different strutjet chamber pressures. The k-e model with and without compressibility corrections, and the SST model 
were used in these CFD calculations. As pressure level of the ejector changes, the flow pattern changes drastically. 
At a chamber pressure of 600 psi, the mixing length was about 0.88 m, from experimental data. The complete 
mixing of the primary and secondary streams is indicated when the side wall pressures recover to the external 
condition, i.e. one ATM. 



Distance, m 


Figure 3. Side-wall pressure comparisons for ejector chamber pressure at 600 psi 


4 

American Institute of Aeronautics and Astronautics 














From these figures, it can be seen that within the mixing zone, the compressibility corrected k-z model shows 
much improved predictions when compared with the standard k-z model and experimental data. However, after the 
complete mixing zone, the wall effect becomes dominant, and the wall function approach used in both k-z models 
(with and without compressibility correction) cannot capture the pressure recovery. On the other hand, the SST 
model performs well in the recovery regions due to its use of k - co formulation at the wall region as seen clearly in 
Figures 3 and 5. 


Three Dimensional Flow Characteristics 

One flow characteristic associated with the non-axisymmetric 
ejector flow field is the three-dimensionality in the strut-base region. 
As depicted in Figure 6, as the primary stream is accelerated through a 
convergent-divergent nozzle to a supersonic Mach 3 jet, the secondary 
stream flows through the strut gap and is then entrained by the 
primary flow. When the secondary flow tries to expand around the 
comer as it passes the strut, a separation region will be formed, 
characterized by low static pressures. The experimental study of 
Reference 1 1 also supplied pressure distribution data at the top wall. A 
full 3-D CFD calculation is required to capture this flow feature. 
Preliminary results using the wall-resolved Baldwin-Lomax 
turbulence model is shown in Figure 7. The good agreement between 
the prediction and data before 0. 1 m and the ability of capturing the 
flow separation by the Baldwin-Lomax model indicated that 
turbulence models capable of resolving near wall regions should be 
used for this 3-D ejector flow field. Compressibility correction to the 
Forsythe et al. 12 , should be further investigated. 


Base Pressure 

TV 
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Secondary Flow 
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Figure 6. Strut-base separation region 
SST model, similar to the approach taken by 



Figure7. 3-D Pressure profiles at centerline of top-wall for ejector chamber pressure 600 psi. 
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IV. SUMMARY 


A RANS numerical modeling study has been carried out to investigate the flow characteristics of the non- 
axisymmetric strut/ejector-based system. In the flow regions dominated by stream mixing, the compressibility 
correction was found to be necessary to capture the decreased spreading rate of compressible flow streams. The wall 
functions used in the k-z model, even with compressibility correction, cannot reproduce wall-dominated phenomena, 
such as pressure recovery. 3D modeling with wall-resolving turbulence models is required to resolve the complex 
strut-base separation features. Full 3D parallel computations, incorporating compressibility corrected SST model 
are underway. These results will form the basis of a future communication. 
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